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Abstract. The M/GI/m/n queueing system with in homogeneous servers 
and the finite number n of waiting spaces is studied. Let A be the customers 
arrival rate, and let be the reciprocal of the expected service time of a 
customer. Under the assumption A = m/i it is proved that the expected 
number of losses during a busy period is the same value for all n > 1 , while in 
the particular case of the Markovian system M /M /m/n the expected number 
of losses during a busy period is for all n > 0. Under the additional 

assumption that the probability distribution function of a service time belongs 
to the class NBU or NWU, the paper establishes simple inequalities for those 
expected numbers of losses in M/GI/m/n queueing systems. 



1. Introduction 

Analysis of loss queueing systems is very important from both the theoretical and 
practical points of view. While the multiserver loss queueing system M /GI /m/Q 
and its network extensions have been intensively studied (see the review paper of 
Kelly [20 , the book of Ross [27^ and references in these sources), the information 
about M/GI/m/n queueing systems (n > 1) is very scanty, because explicit results 
for characteristics of these queueing systems are unknown. (In the present paper, 
for multiserver queueing systems the notation M /GI /ni/n is used, where m denotes 
the number of servers and n denotes the number of waiting places. Another notation 
which is also acceptable in the literature is M/Gl /m/m + n.) 

From the practical point of view, M/GI/m/n queueing systems serve as a model 
for telephone systems, where n is the maximally possible number of calls that can 
wait in the line before their service start. The loss probability is one of the most 
significant performance characteristics. In the present paper, we study the expected 
number of losses during a busy period (the characteristic closely related to the 
stationary loss probability) under the assumption that the arrival rate (A) is equal 
to the maximum service capacity (m^), which seems to be the most interesting 
from the theoretical point of view. 

There are two main reasons for studying this case. 

The first reason is that the case A = m/i is a critical case for queueing systems 
with m identical servers, i.e. the case associated with critically loaded systems. 
The theoretical and practical interest in studying heavily loaded loss systems is 
very high, and there are many results in the literature related to the analysis of 
the loss probability in heavily loaded systems. The asymptotic results for losses in 
heavily loaded single server systems (n — > cxd) such as M/GI/l/n and GI/M/l/n 
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and for associated models of telecommunication systems and dams have been stud- 
ied in [4] , [8] , [9] , [H] , [M] and . Heavy-trafRc analysis of losses in heavily loaded 
niultiserver systems have been provided in |T2], [33], [33] and [35]. The mathemat- 
ical foundation of heavy traffic theory can be found in the textbook of Whitt [32] . 
Although the case A = m/i is idealistic, it enables us to understand the possible 
behaviour of the system in certain cases when the values A and m/z are close and 
approach one another as n increases to infinity. (Obtaining nontrivial results in the 
cases A < m/i and A > m/i is a hard problem, so the analytic investigation of the 
aforementioned asymptotic behaviour as n increases to infinity is difficult.) 

The second reason is that A — m/i is an interesting theoretical case associated 
with an extension of the following non-trivial property of the symmetric random 
walk. Let Xi, X2, Xi, be a sequence of independent and identically 
distributed random variables taking the values ±1 with the equal probability ^. 
Let 5*0 — 0, and S'i+i — Si + ^i+i, i > 0, be a symmetric random walk, and let 
t = T he the first time instant after t = when this random walk returns to zero, 
i.e. St = 0. It is known that the expected number of level-crossings through any 
level n > 1 (or n < — 1) is equal to ^ independently of that level. The mentioning 
of this fact (but in a slightly different formulation) can be found in Szekely [30] . 
and its proof is given in Wolff [37 , p. 411. The reformulation of this fact in terms 
of queueing theory is as follows. Consider M/M/l/n queueing system with equal 
arrival and service rates. For this system, the expected number of losses during a 
busy period is equal to 1 for all n > 0. It has been recently noticed that this property 
holds true for M /GI /1/n queueing systems. Namely, it was shown in several recent 
papers (see Abramov [1], [2], [4], Righter [26|, Wolff [38^), that under mutually equal 
expectations of interarrival and service time, the expected number of losses during 
a busy period is equal to 1 for all n > 0. Further extension of this property to 
queueing systems with batch arrivals have been given in Abramov [5], Wolff |38j 
and Pekoz, Righter and Xia 25 . Applications of the aforementioned property of 
losses can be found in 9 , for analysis of lost messages in telecommunication systems 
and in 11 for optimal control of large dams. Further relevant results associated 
with the properties of losses have been obtained in the paper by Pekoz, Righter and 
Xia [25] . They solved a characterization problem associated with the properties of 
losses in GI /M /1/n queues and established similar properties for M/M/m/n and 
M-^ /M /m/n queueing systems. Recently, a similar property related to consecutive 
losses in busy periods of M /GI /1/n queueing systems has been discussed in [15] . 
It follows from the results obtained in this paper that for M /GI /1/n queueing 
systems with mutually equal expectations of interarrival and service times, the 
expected number of losses in series containing at least fc > 1 consecutive losses 
during a busy period generally depends on n. However, for M/M/l/n queueing 
systems with equal arrival and service rates that expected number of consecutive 
losses during a busy period is the same constant (depending on the value fc) for all 
n > 0. 

The aim of the present paper is further theoretical contribution to this theory of 
losses, now to the theory of multiserver loss queueing systems. On the basis of the 
aforementioned results on losses in M /GI /1/n and M/M/m/n queueing systems 
we address the following open question. Does the result on losses in M/M/m/n 
queueing systems remain true for those M/GI/m/n too'? 
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The answer on this question is not elementary. On one hand, under the as- 
sumption A — mil the expected numbers of losses in M/GI /m/0 and M/GI/m/n 
queueing systems (m > 2 and n > 1) during their busy periods are different. A 
simple example for this confirmation can be built for M/G//2/1 queueing systems 
having the service time distribution G{x) = 1 — pe^^^^ — ge^^^^, p + q = 1. The 
analysis of the stationary characteristics for these systems, resulting in an anal- 
ysis of losses during a busy period, can be provided explicitly. Specifically, the 
structure of the 9x9 Markov chain intensity matrix for the states of the Markov 
chain associated with an M /GI /2/\ queueing system shows a clear difference be- 
tween the structure of the stationary probabilities in M/G//2/1 queues and that in 
M/G//2/0 queues given by the Erlang-Sevastyanov formulae. So, the parameters 
p, q, 111 and /i2 can be chosen such that the expected number of losses during busy 
periods in these two queueing systems will be different. 

On the other hand, the property of losses, which is similar to the aforemen- 
tioned one, indeed holds. The correctness of this similar property for multiserver 
M/GI/m/n queueing systems is proved in the present paper. Namely, we establish 
the following results. 

Let Lm,n denote the number of losses during a busy period of the M/GI/m/n 
queueing system, let A, /i be the arrival rate and, respectively, the reciprocal of the 
expected service time, and let m, n denote the number of servers and, respectively, 
the number of waiting places. We will prove that, under the assumption A — mii, 
the expected number of losses during a busy period of the M/GI/m/n queueing 
system, EL„i^„, is the same for all n > 1, which is not generally the same as that for 
the M/GI/m/0 loss queueing system (when n — 0). In addition, if the probability 
distribution function of the service time belongs to the class NBU (New Better than 
Used), then EL„i_„ = where a constant c > 1 is independent of n > 1. In 

the opposite case of the NWU (New Worse than Used) service time distribution we 
correspondingly have EiLm,n = with a constant c < 1 independent of n > 1 as 
well. (The constant c becomes equal to 1 in the case of exponentially distributed 
service times.) Recall that a probability distribution function S(a;) of a nonnegative 
random variable is said to belong to the class NBU if for all a; > and y > we 
have 'E.{x + y) < S(a;)S(2/), where S(x) = 1 — S(x). If the opposite inequality holds, 
i.e. E{x + y) > S(a;)S(y), then 5(x) is said to belong to the class NWU. 

The proof of the main results of this paper is based on an application of the 
level-crossing approach to the special type stationary processes. The construction 
of the level-crossings approach used in this paper is a substantially extended version 
of that used in the earlier papers by the author (e.g. [1], [3], [6], [7], [10] and [13]) 
and by Pechinkin [23] . It uses modern geometric methods of analysis and involves 
an algebraically close system of processes and a nontrivial construction of deleting 
intervals and merging the ends together with nontrivial applications of the PASTA 
property. 

Throughout the paper, it is assumed that m > 2. (This is not the loss of 
generality since the case m = 1 is known, see [1], [55] and [38).) 

The paper is organized as follows. In Section [2l which is the first part of the 
paper, M/M/m/n queueing systems are studied. The results for M/M/m/n queue- 
ing systems are then used in Section [3l which is the second part of the paper, in 
order to study M/GI/m/n queueing systems. The study in both of Sections [2] 
and [3] is based on the level-crossing approach. The construction of level-crossings 



4 



VYACHESLAV M. ABRAMOV 



for M/M/m/n queueing systems is then developed for M jGl jrajn queueing sys- 
tems as foUows. The stationary processes associated with these queueing systems 
is considered, and the stochastic relations between the times spent in state to — 1 
associated with to — 1 busy servers during a busy period of M/GI /m/n (n > 1) 
and M/GI /ra — 1/0 queueing systems are established. To prove these stochastic 
relations, some ideas from the paper of Pechinkin [24] are involved to adapt and 
develop the level-crossing method for the problems of the present paper. The ob- 
tained stochastic relations are crucial, and they are then used to prove the main 
results of the paper in Section 4. In Section 5, possible development of the results 
for M-^ /GI /ni/n queueing systems with batch arrivals is discussed. 

2. The M/M/m/n queueing system 

In this section, the Markovian M/M/m/n loss queueing system is studied with the 
aid of the level-crossings approach, in order to establish some relevant properties 
of this queueing system. Those properties are then developed for M /GI /m/n 
queueing systems in the following sections. 

Let /(j), 1 < i < n + m -\- \, denote the number of customers arriving during a 
busy period who, upon their arrival, meet j — 1 customers in the system. It is clear 
that /(I) = 1 with probability 1. Let ij.2,- ■ • , tj,f(j) be the instants of arrival 
of these /(j) customers, and let Sj^i, Sj,2,- ■ • , be the instants of the service 

completions when there remain only j — 1 customers in the system. Notice, that 
tn+m+i,k = s«+m+i,fe for aU = 1, 2, . . . , /(n -t- TO -I- 1). 

For 1 < j < + let us consider the intervals 

(2.1) (ij,i, Sj,i], (ij,2, Sj,2], • ■ • , Sj-JO-)]- 
Then, by incrementing index j we have the following intervals 

(2.2) (tj+1,1, Sj+i_i], (<j+i,2, Sj+1,2], • • • , 

Delete the intervals of (|2.2I) from those of (|2.ip and merge the ends, that is each 
point tj+i k with the corresponding point Sj+i k, k — \, 2, /(j + 1) (see Figure 

Then f{j + l) has the following properties. According to the property of the lack 
of memory of the exponential distribution, the residual service time for a service 
completion, after the procedure of deleting the interval and merging the ends as it 
is indicated above, remains exponentially distributed with parameter /imin(j, to). 
Therefore, the number of points generated by merging the ends within the given 
interval (ij,i , Sj.i] coincides in distribution with the number of arrivals of the Poisson 
process with rate A during an exponentially distributed service time with parameter 
/Ltmin(j,TO). Namely, for 1 < j < to — 1 we obtain 

E{/(i + l)|/(j) = l}-^ii / e-^-^j7ie-^''-dx - -. 

Considering now a random number /(j) of intervals (12.11) we have 



(2.3) 



E{/(j + l)|/(j)} = -/(.?). 
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Figure 1 . Level crossings during a busy period in a Markovian system 
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Analogously, denoting the load of the system by p = for m < j < m + n we 
have 

(2.4) mu + = —fU) = pfij)- 

The properties (|2.3p and (|2.4p mean that the stochastic sequence 

(2.5) |/(j + 1) fl min(i, m), J-,+i| , = a{f{l)J{2), . . . , f{j)}, 

forms a martingale. 

It follows from (f231) that for < j < m - 1 

(2-6) E/(, + 1) = 

and for m < j < m -{- n 

(2.7) E/(j + l) = ^^p^-'". 

For example, when p = 1 from (|2.7p we obtain the particular case of the result of 
Pekoz, Righter and Xia |25j: ELm,n=E/(n + m + 1) — for all n > 0, where 
Lm,n denotes the number of losses during a busy period of the M/M/m/n queueing 
system. 

Next, let B{j) be the period of time during a busy cycle of the M/M/m/n 
queueing system when there are exactly j customers in the system. For < j < 
rn + n we have: 



A" 



for < < m — 1, 
.pj-ni^ for m < j < m + n. 



(2.8) AEB(j) = E/(j + 1) 

Now, introduce the following notation. Let Tm,n denote the length of a busy period 
of the M/M/m/n queueing system, let T,„^o denote the length of a busy period of 
the M/M/m/0 queueing system with the same arrival and service rates as in the 
initial M/M/m/n queueing system, and let C,n denote the length of a busy period 
of M/M/\/n queueing system with arrival rate A and service rate [im. From ()2.6p - 
(j2.8p for the expectation of a busy period of the M/M/m/n queueing system we 
have 

n-\-m m — 1 — i \m— 1 " 

(2-9) Er^,„ J2 E^C?) = E -TT + -nr^ E p"- 

j=l j = i ^ j=0 

In turn, for the expectation of a busy period of the M/M/m/0 queueing system we 
have 



(2.10) Er„,o = EE^W = 



where (|2.10p is the particular case of (|2.9p where n = 0. 

It is clear that Tm,n contains one busy period Tm-i,o, where the subscript m — 1 
underlines that there are m — 1 servers, and a random number of independent busy 
periods, which will be called orbital busy periods. Denote an orbital busy period by 
C,n- (It is assumed that an orbital busy period starts at instant when an arriving 
customer finds m—1 servers busy and occupies the mth server, and it finishes 
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at the instant when after a service completion there at the first time remain only 
TO — 1 busy servers.) Therefore, denoting the independent sequence of identically 
distributed orbital busy periods by c},i \ C,n\---, we have 



(2.11) 2^jn,n— Tm-1,0 + C 

1=1 

where k is the random number of the aforementioned orbital busy periods and = 
means an equality in distribution. It follows from ()2.9p . (|2.10|) and ()2.11|) 

(2-12) E^CW^^E/'^- 

i=i ^ j=0 

On the other hand, the expectation of an orbital busy period is 

1 " 

(this can be easily checked, for example, by the level-crossings method [I], [B] and 
an application of Wald's identity :16 , p. 384), and we obtain 

^m — 1 

(2.13) Ek=- -. 

^ ' (to - l)!/x™-i 

Thus, Ek coincides with the expectation of the number of losses during a busy 
period in the M/M/m/0 queueing system. In the case p = 1 we have Ek = ^L^. 

3. M/GI/m/n queueing systems 

In this section, the inequalities between the times spent in the state to — 1 in the 
M/GI/m/n (n > 1) and AI/GI /m/0 queueing systems during their busy periods 
are derived. 

Consider two queueing systems: M/GI/m/n [n > 1) and M/GI/m/0 both 
having the same arrival rate A and probability distribution function of a service 
time G{x), i — xdG{x) < oo. Let T„ „(to — 1) denote the time spent in the 
state TO — 1 during its busy period (i.e. the total time during a busy period when 
TO— 1 servers are occupied) of the M/GJ/m/n queueing system, and let Tm_o("^ — 1) 
have the same meaning for the M/GI/m/0 queueing system. 

We prove the following lemma. 

Lemma 3.1. Under the assumption that the service time distribution G(x) belongs 
to the class NBU (NWU), 

(3.1) T„i,n{m - 1) >sf {resp. <st) T,nfi{m - 1), 

Proof. . The proof of the lemma is relatively long. In order to make it transparent 
and easily readable we strongly indicate the steps of this proof given by several 
propositions (properties). There are also six figures (Figures 2-7) illustrating the 
constructions in the proof. Each of these figures contain two graphs. The first 
(upper) of them indicates the initial (or intermediate) possible path of the process 
(sometimes two-dimensional) , while the second (lower) one indicates the part of the 
path of one or two-dimensional process after a time scaling or specific transforma- 
tion (e.g. in Figure 5). Arc braces in the graphs indicate the intervals that should 
be deleted and their ends merged. 
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Two-dimensional processes are shown as parallel graphs. For example, there are 
two parallel processes in Figure 3 which are shown in the upper graph, and there 
are two parallel processes which are shown in the lower graph. The same is in 
Figures 4, 6 and 7. 

For the purpose of the present paper we use strictly stationary processes of order 
1 or strictly 1- stationary processes. Recall the definition of a strictly stationary 
process of order n (see [53], p. 206). 

Definition 3.2. The process ^(t) is said to be strictly stationary of order n or 
strictly n-stationary, if for a given positive n < oo, any h and ti, t2,- ■ ■ , tn the 
random vectors 



have identical joint distributions. 

If n = 1 then we have strictly 1-stationary processes satisfying the property: 



The probability distribution function P{^(i) < x} in this case will be called limiting 
stationary distribution. 

The class of strictly 1-stationary processes is wider than the class of strictly 
stationary processes, where it is required that for all finite dimensional distributions 



pmti), at2), . . . , e(ife)) e = pmh + h),at2 + h),..., atk + h}} e s,}. 



for any h and any Borel set Bk C M.^. The reason of using strictly 1-stationary 
processes rather than strictly stationary processes themselves is that, the operation 
of deleting intervals and merging the ends is algebraically close with respect to 
strictly 1-stationary processes, and it is not closed with respect to strictly stationary 
processes. The last means that if ^(t) is a strictly 1-stationary process, then for 
any h > and arbitrary to the new process 



is also strictly 1-stationary and has the same one-dimensional distribution as the 
original process £,{t). The similar property is not longer valid for strictly stationary 
processes. If ^{t) is a strictly stationary process, then generally (,i{t) is not strictly 
stationary. 

In the following the prefix 'strictly' will be omitted, so strictly stationary and 
strictly 1-stationary processes will be correspondingly called stationary and 1- 
stationary processes. 

Let us introduce m independent and identically distributed stationary renewal 
processes (denoted below Xm(t)) with a renewal period having the probability dis- 
tribution function G{x). 

On the basis of these renewal processes we build the stationary m-dimensional 
Markov process Xm(t) = {^i{t), ^2{t), ■ . ■ , £,m{t)}, the coordinates ^k{t), k =1,2,. . . , 
m of which are the residual times to the next renewal times in time moment t, 
following in an ascending order. 

Let us now consider the two (m -f l)-dimensional Markov processes correspond- 
ing to the M/GI/m/n (n > 1) and M/GI/m/0 queueing systems, which are de- 
noted by ym,n(t) and Ymfiit). Let Qm,n{t) denote the stationary queue-length 




P{^{t)<x} = P{^{t + h)<x}. 



^{t + h), ift>to 
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process (the number of customers in the system) of the M/GI/m/n queueing sys- 
tem, and let Qm,o(t) denote the stationary queue- length process corresponding to 
the M/GI/m/Q queueing system. We have: 

Ym.nit) = {r7^'"^(t),ryf'"^(t),...,^L"'")(i),Q™,„(t)}, 

where 

\^m-iy„,„(t) + ll''^ '/m-i/,„,„(t)+2'^''''' ' • ' ' '/m y'^) f 

are the ordered residual service times corresponding to i'm,n{t) = min{r7i, Qm,n{t)} 
customers in service in time t, and 

{^("-) {tl ryf-") (t) , . . . , (0 } = {0, 0, . . . , 0} 

all are zeros. Analogously, 

y™,o(t) = {77^'°)(t),r7f'°^(t),...,77^-°)(i),g,„,o(t)}, 

only replacing the index n with 0. 

Let us delete all time intervals of the process ym,n{t) related to the M/GI/m/n 
queueing system {n > 1) where there are more than m — 1 or less than to — 1 
customers and merge the ends. Remove the last component of the obtained process 
which is trivially equal to to — 1. Then we get the new (to — l)-dimensional Markov 
process (in the following the prefix 'Markov' will be omitted and only used in the 
places where it is meaningful): 

y™-i.„(t) = {^^"^(^),r7^")(^),...,57i^J?(^)}, 

the components of which are now denoted by hat. All of the components of this 
vector are 1-stationary, which is a consequence of the existence of the limiting 
stationary probabilities of the processes r]'™''"^\t), j = 1, 2, . . . , to (e.g. Takacs [3T| ) 
and consequently those of the processes rf™''"'^ (i), j = 1, 2, . . . , to. The joint limiting 
stationary distribution of the process ym-i.n{t) can be obtained by conditioning of 
that of the processes ym,n(t) given Qm,nit) — m ~ I. 

The similar operation of deleting intervals and merging the ends, where there 
are less than m — 1 customers in the system, for the process ym-i.o(i) is used. We 
correspondingly have 

We establish the following elementary property related to the M/GI/l/n queue- 
ing systems, n=0,l,. . . 

Property 3.3. 

(3.2) P {v^^'"'\t) G Bi I gi,„(i) > l} = P{xi(<) e Si}, 

for any Borel set Bi C . 

Proof. Delete all of the intervals where the server is free and merge the correspond- 
ing ends (see Figure 2). Then in the new time scale, the processes all are structured 
as a stationary renewal process with the length of a period having the probability 
distribution function G{x). Therefore (|3.2p follows. □ 
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Residual service times in tlie original M/GI/l/n queueing system 




Figure 2. Residual service times for the original and scaled pro- 
cesses of the M/GI /1/n queueing system. 
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In order to establish similar properties in the case m = 2 let us first study the 
properties of 1-stationary processes and explain the construction of tagged server 
station which is substantially used in our construction throughout the paper. 

Properties of 1-stationary processes. Recall (see Definition 13.21) that if S^it) is a 
1-stationary process, then for any h and tg the probability distributions of ^(to) 
and ^(^0 + h) are the same. The result remains correct (due to the total proba- 
bility formula) if h is replaced by random variable •& with some given probability 
distribution, which is assumed to be independent of the process ^. Namely, we 
have: 



(3.3) 



/OO 
P{e(io + h)< x]AV{'d < h} 
-OO 

/OO 
dp{d < h} 
-oo 



= PU(io) <4- 
That is, ^(^o) and ^(^o + have the same distribution. 

The above property will be used for the following construction of the sequence 
of 1-stationary processes ^(i)(t), . . . ^ having identical one-dimensional dis- 

tributions. 

Let ^^°Ht) = £,{t) be a 1-stationary process, let ti be an arbitrary point, and 
let t?i be a random variable with some given probability distribution, which is 
independent of the process C^^HO- Let us build a new process 

C'^H^) as follows. 

Put 



(3.4) e^it) 



\^^°^{t + di), foraUt>ti. 



Since the probability distributions of ^{t) and + ^i) are the same for all t > ti, 
then the processes £_{t) and C'^H^) have the same one-dimensional distributions, 
and '^'■^-'(i) is a 1-stationary process as well. 

With a new point t2 and a random variable i?2 , which is assumed to be indepen- 
dent of the process (t) and random variable (therefore, it is also independent 
of the process by the same manner one can build the new 1-stationary 

process 

e^2)(i). Specifically, 



(3.5) e'\t) 



£,^^Ht), foraUt<t2, 
+ foraUt>i2. 



The new process 

^(2)(i) has the same one-dimensional distribution as the processes 
^^°\t) and C^^H^)- The procedure can be infinitely continued, and one can ob- 
tain the infinite family of 1-stationary processes, having the same one-dimensional 
distribution. 

The points ti, t2,...i'n the above construction are assumed to be some fixed 
(non-random) points. However, the construction also remains correct in the case 
of random points to, ti,. . . of Poisson process, since according to the PASTA prop- 
erty [36] the limiting stationary distribution of a 1-stationary process in a point 
of a Poisson arrival coincides with the limiting stationary distribution of the same 
1-stationary process in an arbitrary non-random point. Furthermore, the aforemen- 
tioned property of process remains correct when the random points to,ti,... are the 
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points of the process which is not necessarily Poisson but belongs to the special class 
of processes that contains Poisson. In this case the property is called ASTA (e.g. 

l-stationary Poisson process. Consider an important particular case when the 
process C(t) is Poisson. Let ^(^'(i) = £_{t). Then the process ^(i)(i) that obtained 
by (|3.4p is no longer Poisson. Its limiting stationary distribution is the same as 
that of the original process ^ (t) , but the joint distributions of this process given in 
different points s and t distinguish from those of the original process ^(t). 

The process ^'^'(t) will be called l-stationary Poisson process or simply 1- 
Poisson. Clearly, that the further processes such as S}-"^^ (t), ^^^^ (t), ... that obtained 
similarly to the procedure in p.4p . (|3.5|) all are 1-Poisson with the same limiting 
stationary distribution. According to the above construction, a 1-Poisson process is 
obtained by deleting intervals and merging the ends of an original Poisson process. 
Therefore, a sequence of 1-Poisson arrival time instants is a scaled subsequence of 
those instants of the ordinary Poisson arrivals. Hence, for 1-Poisson process the 
ASTA property is satisfied, i.e. 1-Poisson arrivals see time averages exactly as those 
Poisson arrivals. 

Tagged server station. Consider a stationary queueing system M/GI/m/n, 
which is referred to as main server station, and in addition to this queueing system 
introduce another one containing a server station in order to register specific ar- 
rivals, for example losses or, say, customers waiting their service in the main system. 
This server station is called tagged server .station. The main idea of introducing 
tagged server stations is to decompose the main system as follows. Assume that 
along with a Poisson stream of arrivals of customers occupying servers in the main 
system, there is another stream of arrivals of customers in the tagged server system. 
For instance, the losses in the main system can be supposed to occupy the tagged 
server station. Although the stream of these losses is not Poisson (see e.g. [211, p. 
83 or [20], p. 320), it is shown later that it is 1-Poisson. Therefore, the original 
system is decomposed into smaller systems with the same (1-Poisson) type of input 
stream. It is worth noting that only one dimensional distributions of 1-Poisson 
process are the same for all of them that generated similarly to the procedure in 
p.4p . (|3.5p . However, the two-dimensional distributions are distinct in general. 

In fact, applications of a tagged server station is wider than that, and its aim 
is a proper decomposition of the original system into the main and tagged systems 
for further study of the properties of losses. 

Another idea of using tagged server stations is a proper application of the ASTA 
property as follows. At the moment of arrival of a customer in the tagged server 
station, the stationary characteristics in the main server station remain the same. 
Specifically, the distributions of residual service times in servers of the main station 
at the moment of arrival of a customer in the tagged station coincides with the 
usual stationary distributions of these residual service times. 

Let us now formulate and prove a property similar to Propertv 13.31 for m = 2. 
We have the following. 

Property 3.4. For the M/GI/2/0 queueing system we have: 



(3.6) 



P{yi,o(t) e Bi} = P{xi(t) e B,}, 
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Proof. In order to simplify the explanation in this case, let us consider two auxiliary 
stationary one-dimensional processes Ci,o(^) and C2,o(^)- The first process describes 
a residual service time in the first server, and the second one describes a residual 
service time in the second server. If the ith server [i = 1, 2) is free in time i, then 
we set Ci.o(^) ■— 0. 

Our further convention is that the first server is a tagged server. We assume 
that if at the moment of arrival of a customer both of the servers are free, he/she 
occupies the first server. Clearly that this assumption is not a loss of generality. For 
instance, if we assume that both of the servers are equivalent and can be occupied 
with the equal probability i, then an occupied server (let it be the first) can be 
called tagged. In another busy period start an arriving customer can occupy the 
second server. It this case, nothing is changed if the servers will be renumbered, 
and the occupied server will be numbered as first and called tagged. 

Our main idea is a decomposition of the stationary M /GI /2/Q queueing system 
into two systems and study the properties of stationary (1-stationary) processes 
Ci,o(0 and C2,o(i)- The arrival stream to the tagged system is Poisson, so the first 
system is M/G//1/0, while the second one is denoted •/G//1/0, where • in the 
first place of the notation stands for the input process in the second system, which 
is the output (loss) stream in the first one. Clearly, that an arriving customer is 
arranged to the second queueing system if and only if at the moment of his/her 
arrival the tagged system is occupied. Therefore, let us delete all the intervals when 
the tagged system is empty and merge the ends. In this case, the tagged system 
becomes an ordinary renewal process, and the stream of arrivals to the second 
queueing system becomes 1-Poisson rather then Poisson (because after deleting 
intervals and merging the ends in the new time scale the original Poisson process is 
transformed into 1-Poisson). Therefore the second system now can be re-denoted 
by M/G// 1/0, where M in the first place of the notation stands for 1-Poisson input 
and replaces the initially written symbol •. 

Thus, the M/G//2/0 queueing system is decomposed into the Af/G//l/0 and 
M/G//1/0 queueing systems. Clearly, that without loss of generality one can 
assume that the original arrival stream is 1-Poisson rather than Poisson, i.e. the 
original queueing system is M/G//2/0, and it is decomposed into two M/G// 1/0 
queueing systems. The last note is important for the further extension of the result 
for the systems M /GI /m/Q (or generally M /GI /m/Q) having m > 2 servers. 

Let T be the time moment when an arriving customer occupies the tagged server 
station. According to the ASTA property, 

(3.7) P{C2,o(t) < x} = P{C2,o(i) < x}, 

where t is an arbitrary fixed point, and the probability distribution function of 
C2,o(0 in this point coincides with the distribution of residual service time in specific 
M I GI /I/O system with some specific value of parameter of 1-Poisson process, which 
is not important here. On the other hand, the process C2,o(0 is stationary and 
Markov. Therefore from p. 71) for any /i > we have 

(3.8) P{C2,o(r + h)<x]= P{C2,o(i + h) < x] = P{C2,o(i) < x]. 

Let X denotes the service time of the customer, who arrives at the time moment 
r occupying the tagged server station. Our challenge is to prove that 



(3.9) P{C2.o(r + x)<x}^ P{C2,o(t) < a:|Ci,o(i) > 0}. 
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Instead of the original processes C,i.o{t), i — 1,2, consider another processes Ci,o(i), 
which are obtained by deleting the intervals where the tagged server is free, and 
merging the ends. Then, Ci,o(i) is a renewal process, and the 1-stationary process 
C2,o(t) and the random variable x (the length of a service time in the tagged server 
that starts at moment r) are independent. Hence, for any event {% = h} according 
to the properties of 1-stationary processes we have 

(3.10) P{C2,o(t + x) < = h}^ P{C2,o(r) < x}, 
and, due to the total probability formula from (|3.10l) we have 

(3.11) P{C2,o(t + x)<x} = P{C2,o(t) < x}. 

The only difference between p. lip and the basic property p.3p is that the time 
moment r is random, while to is not. However keeping in mind p.Sp . this modified 
equation p. lip follows by the same derivation as in p.3p . 
Hence, from p. lip . 

(3.12) P{C2,o(t + x)<^}^ P{C2,o(t) <x}^ P{C2,o(i) < x\CMt) > 0}, 

and since P{C2,o(''" + x) ^ ^} — P{C2,o(''' + x) ^ ^} ^e finally arrive at p.Op . 
As well, noticing that 

P{C2,o(r) < x} = P{C2,o(t) < x}, 
from p. lip and p.7p we also have 

(3.13) P{C2,o(t + x)<x}^ P{C2,o(t) <x}= P{C2,o(i) < x}. 
Similarly to p.l3p one can prove 

(3.14) P{Ci,o(t + x)<x}^ P{Ci,o(t) < x} ^ P{Ci,o(t) < x}, 

where r is the moment of arrival of a customer, who at this moment r occupies the 
second server, and x denotes his/her service time. Relations (I3.14p can be proved 
with the aid of the same construction of deleting intervals and merging the ends, 
but now in the second server. So, combining p.l3p and p.l4p we arrive at the 
following fact. In any arrival or service completion time instant in one server, the 
residual service time in another server has the same stationary distribution. 
This fact is used in the constructions below. 

Now consider the stationary M/GI/2/0 queueing system, in which both servers 
are equivalent in the sense that if at the moment of arrival of a customer both servers 
are free, then a customer can occupy each of servers with the equal probability ^. 
In this case the both of the processes Ci,o(0 C2,o(0 have the same distribution. 

Let us delete the time intervals where the both servers are simultaneously free, 
and merge the corresponding ends (see Figure 3). The new processes are denoted by 
Ci,o{t) and C2,o(i), and both of them have the same equivalent distribution. (We use 
the same notation as in the construction above believing that it is not confusing 
for readers.) This two-dimensional 1-stationary process characterizes the system 
where in any time t at least one of two servers is busy. Consider the event of arrival 
of a customer in a stationary system at the moment when only one of two servers 
is busy. Let r* be the moment of this arrival, and let r** denote the moment of the 
first service completion in one of two servers following after the moment r* . Then, 
at the endpoint t** of the interval [t*,t**) the distribution of the residual service 
time will be the same as that at the moment r* (due to the established fact that at 
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the end of a service completion in one server, the distribution of a residual service 
time in another server must coincide with the stationary distribution of a residual 
service time and due to the fact that both servers are equivalent.) 

The additional details here are as follows. There can be different events asso- 
ciated with the points r* and t**. For example, at the moment r* an arriving 
customer can be accepted by one of the servers, while the service completion at the 
moment t** can be either in the same server of in another server. If time moments 
T* and T** are associated with the same server (for example, the moment of service 
start and service completion in the first server) then we speak about residual service 
times in another server (in this example - the second server). If time moments r* 
and T** are associated with different servers (say, r* is the service start in the first 
server, but r** is the service completion in the second one), then we speak about 
residual service times in different servers (in this specific case we speak about resid- 
ual service time in the second server at the time moment r* and residual service 
time in the first server at the time moment r**). However, according to the earlier 
result, it does not matter which specific event of these mentioned occurs. The only 
fact, that the stationary distribution of a residual service time in a given server 
must be the same for all time moments of arrival and service completion occurring 
in another server and vice versa, is used. 

Deleting the interval [t*, t**) and merging the ends r* and r** (see Figure 4) 
we obtain the following structure of the 1-stationary process yi,o(^)- 

In the points where idle intervals are deleted and the ends are merged we have 
renewal points: one of periods is finished and another is started. In the other 
points where the intervals of type [t* , t**) are deleted and their ends are merged 
we have the points of 'interrupted' renewal processes. In this 'interrupted' renewal 
process the point r* is a point of 1-Poisson arrival, and, according to ASTA, the 
distribution in this point in the server that continue to serve a customer coincides 
with the stationary distribution of the residual service time. In the other point 
T**, which is the point of a service completion, the distribution in this point in the 
server that continue to serve a customer coincides with the stationary distribution 
of a residual service time as well. Therefore, in the point of the interruption (which 
is a point of discontinuity) the residual service time distribution coincides with 
the stationary distribution of a residual service time, i.e. with the distribution of 
Xi(i). (Notice, that the intervals of type [t* ,t**) are an analogue of the intervals 
[si.fej^i.fe) considered in the Markovian case in Section[5J) 

By amalgamating the residual service times of the first and second servers given 
in the lower graph in Figure 4, one can built a typical one-dimensional 1-stationary 
process yi.o(i), the limiting stationary distribution of which coincides with that of 
Xi(t). (see Figure 5). 

Therefore the processes yi.o(i) and Xi(t) have the identical one-dimensional dis- 
tribution, and relation (I3.6|) follows. □ 

Let us develop Propertv l3.4l to the case m — 3 and then to the case of an arbitrary 
m > 1 for the M/GI/m/0 queueing systems. Namely, we have the following. 

Property 3.5. For the M/GI/m/0 queueing system we have: 



(3.15) 



P{ym-l,o(0 e B„^-l} = P{x„._i(t) G B„^-l}, 



where Bm-i is an arbitrary Borel set ofW' 



m — 1 
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Residual service times in tlie original M/GI/2/0 queueing system 




Residual service times of the scaled process in the M/GI/2/0 queueing system 
after deleting the intervals when both of the servers are idle, and merging the 
ends 




Figure 3. Residual service times for the original and scaled pro- 
cesses of the M/GI/2/Q queueing system after deleting intervals 
where both of the servers are free, and merging the ends. 
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Residual service times of tlie scaled process in tlie M/GI/2/0 queueing system 
after deleting the intervals when both of the servers are idle, and merging the 
ends 




Residual service times of the scaled process in the M/GI/2/0 queueing system 
after deleting the intervals vchen both of the servers are idle and both are busy, 
and merging the ends 




Figure 4. Residual service times for the original and scaled pro- 
cesses of the M/G//2/0 queueing system after deleting intervals 
where both of the servers are free and both are busy, and merging 
the ends. 
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Figure 5. A typical 1-stationary process of residual service times 
is obtained by amalgamating residual service times of the first and 
second servers, rj* and t| are the points where the intervals of type 
[t*,t**) are deleted, and the ends are merged 
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Proof. The proof will be concentrated in the case to = 3 for the 1-stationary process 
y2,o(^), which is associated with the paths of the M/GI/3/0 queueing system where 
only two servers are busy. Then the result will be concluded for an arbitrary m > 2 
by induction. 

Prior studying this case, we first study the specific case of the M/GI/2/0 queue- 
ing system by considering the paths when the both servers are busy. Then using 
the arguments of the proof of Propertv l3. 41 enables us to extend that specific result 
related to the M/GI/2/0 queueing systems to the 1-stationary process y2,o(i) of 
the M/GI/3/0 queueing system. 

As in the proof of Propertv 13.41 in the specific case of the M/GI /2/Q queueing 
system considered here, we will study the stationary one-dimensional processes 
Ci,o(i) and C2,o(i)- However the idea of the present proof generally differs from 
that of the proof of Property 13.41 Here we do not call the first (or second) server 
a tagged server station to use decomposition. We simply use the fact established 
in the proof of Property 13.41 that at the moment of arrival or service completion 
of a customer in one server, the distribution of a residual service time in another 
server will coincide with the stationary distribution of a residual service time in 
this server. (The same idea has been used in the proof of Propertv 13.41 ) 

The present proof explicitly uses the fact that the class of 1-stationary processes 
is algebraically closed with respect to the operations of deleting intervals and merg- 
ing the ends, which was mentioned before. 

Let us delete the idle intervals of the process Ci,o(^) and merge the ends. Then 
we get a stationary renewal process as in the above case m — 1 (Propertv 13. 3p . 

After deleting the same time intervals in the second stationary process C2,o(^) and 
merging the ends, the process will be transformed as follows. Let t* be a moment of 
1-Poisson arrival when a customer occupies the first server. (Recall that owing to 
the known properties of 1-Poisson process, the stream of arrival to each of i servers 
(i = 1,2) is 1-Poisson.) Then, according to the ASTA property, C2.o(i*) = C2,o(i) 
in distribution. Therefore after deleting all of the idle intervals of the second server 
and merging the ends, after the first time scaling (i.e. removing corresponding time 
intervals, see Figure 6) instead of the initial 1-stationary process C2,o(^) we obtain 
the new 1-stationary process with the equivalent one-dimensional distribution. This 
process is denoted by C2,o(i)- 

Notice, that the process C2,o(^) is obtained from the process C2,o(i) by construct- 
ing a sequence of 1-stationary processes described above. 

Then we have the two-dimensional process the first component of which is xi (t) 
and the second one is C,2,o{t)- For our convenience this first component is provided 
with upper index, and the two-dimensional vector looks now as |x^^^(i), C2,o(0|- 

Let us repeat the above procedure, deleting the remaining idle intervals of the 
second server and merging the ends. We get the 1-stationary process being equiv- 
alent in the distribution to the stationary renewal process Xi(i), which is denoted 
now xj^-* (t) . 

Upon this (final) time scaling the first process x^^'' (t) is transformed as follows. 
Let t** be a random point of 1-Poisson arrival when the second server is occupied. 
Applying the ASTA property once again, for the first component of the process we 
obtain that T6^\t**) coincides in one-dimensional distribution with yL^^\t). Thus, 
after deleting the entire idle intervals and merging the ends, we finally obtain the 
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The processes ^i o(Oand (^2,0 W 




The transformed processes t^, „ (t) and t^^ after deleting the intervals and 
merging the ends 



▲ 






\ 




► 

► 



Figure 6. The dynamic of time scaling for a queueing system 
with two servers after deleting the idle intervals in the first server 
and merging the ends 
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two-dimensional process {t),yi\ [t) | each component of which has the same 

one-dimensional distribution as this of the process Xi (i) . The dynamic of this time 
scaling is shown in Figure 7. 

For our further purpose, the independence of the processes xj^'' (i) and x^^"* (t) 
is needed. The constructions in this paper enables us to prove this independence. 
However, the independence of 'x!"i\t) and 'x!"i\t) follows automatically from the 
known results by Takacs [31 and the easiest way is to follow a result of that 
paper. Namely, it follows from formulae (6) and (7) on page 72, that the joint 
conditional stationary distribution of residual service times given that k servers 
are busy coincides with the stationary distribution of Xk{t), which in turn is the 
product of the stationary distributions of Xi(t). In particular, 

(3.16) P [x.\'\t) < XI, xf)(0 <a;2} =P{xi(0 <.Ti}P{xi(i) <a;2}. 

Now, using the arguments of the proof of Propertv 13.41 one can easily extend the 
result obtained now for M/GI/2/Q queueing system to the M/GI/'i/Q queueing 
system, and thus prove (I3.15P for the M/GI/'i/G queueing system. 

Similarly to the proof of Propertv 13.41 let us introduce the processes Ci,o(i), 
(2,0(0 ^i^d C3,o(0 of the residual service times in the first, second and third servers 
correspondingly. These processes all are assumed to have the same stationary 
distribution of residual times, which respects to the scheme where an arriving cus- 
tomer occupies one of available free servers with equal probability. Let us call a 
server of the M/G//3/0 queueing system that occupied at the moment of busy 
period start a tagged server station. So, we decompose the original system into the 
M /GI/2/Q and tagged queueing system M/GI/l/O. However, it is shown above 
that M /GI/2/Q can be decomposed into two M /G//1/0 queueing systems, where 
after the procedure of deleting idle intervals and merging the ends we obtain the 
process having the same stationary distribution as that of the process X2(<). This 
stationary distribution remains the same in all random points of arrivals and service 
completions in the tagged service station. So, after deleting intervals and merging 
the ends in the tagged service station, in a new time scaling we arrive at the same 
stationary distribution as that of the process X3(t). So, the result for m = 3 follows. 

This induction becomes clear for an arbitrary m > 2 as well, where the original 
M /GI /m/Q system can be decomposed into the M /GI /m — 1/0 queueing system 
and a tagged server station M /GI □ 

Now we will establish a connection between the processes ym-i,n(i), ym-i,o(0 
and Xm_i(t). We start from the case m = 2. 

Property 3.6. Under the assumption that the probability distribution function 
G{x) belongs to the class NBU (NWU) we have 

(3.17) P{yi,n(0 <x}< {resp. >) P{xi(i) < x}. 

Proof. Along with the 1-stationary processes yi,n(i), let us introduce another 1- 
stationary processes Y2,„(t). This last process is related to the same M /GI /2/n 
queueing system as the process yi,n(t), and is obtained by deleting time intervals 
when there are more than two or less than two customers in the system, and merging 
the ends. 
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Figure 7. The dynamic of time scaling for a queueing system with 
two servers after deleting the idle intervals in the second server and 
merging the ends 
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Using the same arguments as in the proof of Propertv l3.5l one can prove that the 
components of this process are generated by independent 1-stationary processes 
and having the same distribution. Indeed, involving as earher in the proof of 
Property 13.41 the processes Ci {t) and C2 [t) having the same distribution, one can 
delete intervals where the system is empty and merge the ends. Apparently, the 
new processes Ci (t) and ^2 {i) obtained after this procedure have the same stationary 
distribution. (However, it is shown later that the limiting stationary distribution 
of these one-dimensional processes differs from such the distribution obtained after 
the similar procedure for the M/G//2/0 system, and, therefore, its one-dimensional 
distribution distinguishes from that of the process xi(t).) 

Let us go back to the initial process y2.n{t), n > 1, to delete the time intervals 
where the both servers are free and merge the corresponding ends. We also remove 
the last component corresponding to the queue-length Q2.n{t)- (The exact value of 
the queue-length Q2.n{t) is irrelevant here and is not used in our analysis.) In the 
new time scale we obtain the two-dimensional process y2,n{i)- 

Similarly to the proof of relation p.6|) we have time moments r* and r**. The 
first of them is a moment of arrival of a customer at the system with one busy 
and one free server, and the second one is the following after r* moment of service 
completion of a customer when there remain one busy server only. The time interval 
[r*, T**) is an orbital busy period. (The concept of orbital busy period is defined in 
Section [5] for Markovian systems. For M/GI/m/n queueing systems this concept 
is the same.) It can contain queueing customers waiting for their service. Let 
^begin ^ moment of arrival of a customer during the orbital busy period [r*, r**) 
who occupies a waiting place, and let be the following after i^ogm jj^oment of 
time when after the service completion the queue space becomes empty again. A 
period of time t'^'^'^) is called queueing period. (Note that for the M/GI/m/n 

queueing system, the intervals of type [r*, r**) are an analogue of the intervals 
of type [sm-i,k,im-i,k) in the Markovian queueing system Ad /M /m/n, and the 
intervals of type t'^"'^) are an analogue of the intervals of type [sm.fc, ^m,fc) in 

that Markovian queueing system M /M /m/n.) 

All customers of queueing periods^ i.e. those arrived during orbital busy period 
can be considered as customers arriving in a tagged server station. At the moment 
of t^°s'"j which is an instant of a Poisson arrival, the two-dimensional distribution 
of the random vector y2,n(i'^''^'") coincides with the stationary distribution of the 
random vector Y2,„(i). However, in the point t"'^'^, the probability distribution 
of y2,n{t°'^'^) is different from the stationary distribution of Y2,„(i), because this 
specific time instant coincides with a service beginning in one of servers of 
the main system. Therefore, deleting the interval [£°°^™ ^t°'^'^) and merging the end 
leads to the change of the distribution. 

More specifically, at the time instant t°'^'^ one of the components of the vector 
Y2,n(i), say the first one, is a random variable having the probability distribution 
G{x). Then, another component, i.e. the second one, because of the aforemen- 
tioned properties of 1-stationary processes, coincides in distribution with Ci^„ (or 
C2,n), which is a component of the stationary process Y2,„(t). Indeed, let customers 
arriving in a busy system and waiting in the queue be assigned to the tagged server 
station. At the moment of 1-Poisson arrival ^bogm ^ customer in the tagged server 
station, the two-dimensional Markov process associated with the main queueing 
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system has the same distribution as the vector Y2,n(t), i.e. two-dimensional distri- 
bution coinciding with the joint distribution of and C2.n- Then, at the moment 
of the service completion t"^'^, which coincides with the moment of service com- 
pletion in one of two servers, the probability distribution function of the residual 
service time in another server, where the service is being continued, coincides with 
the distribution of a component of the vector Y2,,i(t), i.e. with the distribution of 

If the probability distribution function G{x) belongs to the class NBU, then the 
1-stationary process y2,n{t) satisfies the property y2,n{t^'^^^'^) <st y2,n(^''"^). If 
G{x) belongs to the class NWU, then the opposite inequality holds: y2,Ti(i°"'^) <st 
y2,n(i'^°^'")- (The stochastic inequality between vectors means the stochastic in- 
equality between their corresponding components.) 

The above stochastic inequalities are between random values of the process 
y2,n{t) in the different time instants t^^s^a g^^^j^ ^cnd^ q^^. f^j-^^i^g]. ^^gj^ jg compare 
two different processes y2,n(i) and y2,o(0- The first of these processes is associated 
with the M/GI/m/n queueing system, while the second one is associated with the 
M/GI /m/0 queueing system. The idea of comparison is very simple. Suppose 
that both queueing system are started at zero, i.e. consider the paths of these 
system when the both of them are not in steady state, and compare the Markov 
processes associated with these system. For the non-stationary processes we will 
use the same notation y2,n 

(t) and y2,o{t) understanding that it is spoken about 
usual (not stationary) Markov processes. The notation for time moments such as 
^begin g^j-^^ ^ond associatcd with these usual (i.e. non-stationary) processes 

as well. We will consider the Markov processes associated with M/GI/2/n and 
M/GI /2/0 queueing systems on the same probability space. In the time interval 
|-Q^^begm^ the paths of the Markov processes y2,n(i) and y2,o(0 coincide (n ^ 0). 
However, after deleting the interval [t^°s^^^ , t°'^'^) and merging the ends, then in the 
end point t'='<'s™ the values of the processes y2,n{t) and y2,o(i) will be different. 
Indeed, in the case of the process y2,o{t), which is associated with the M/GI/m/0 
queueing system, ^bcgm .^^^ ^cnd same point, and the value of Markov pro- 

cesses will be the same after replacing the points i^cgm However, in the 

case of the process y2,n{t) associated with M/GI/m/n queueing system, the values 
in these points will be different with probability 1, and consequently, because of 
the inequality y2,„(i'''=s'") <.t y2,n(t™^) we have y2.o(i'"^*5'") <st y2,nit''"^'") in the 
case when G{x) belongs to the class NBU. If G{x) belongs to the class NWU, we 
have the opposite inequality: y2,n(i''°^''') <st y2,o(i''''^'")- 

Therefore, after deleting all the intervals of the type [t^°^^'^, t^"'*) from the orig- 
inal Markov process we obtain new Markov process, and in the case when G{x) 
belongs either to the class NBU or to the class NWU one can apply the theorem of 
Kalmykov [17] (see also ,19 J to compare these two Markov processes. In the case 
where G{x) belongs to the class NBU, all the path of the Markov process, associated 
with M/GI /2/n is not smaller (in stochastic sense) than that path of the Markov 
process, associated with M/GI/2/0. If G{x) belongs to the class NWU, then the 
opposite stochastic inequality holds between two different Markov processes. Ap- 
parently, the same stochastic inequalities remain correct if we speak about station- 
ary Markov processes. Nothing is changed if we let t to increase to infinity and 
arrive at stationary distributions. So, under the assumption that G{x) belongs to 
the class NBU, for the stationary processes we obtain yi,o(0 ^st yi,n{t)- In other 
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words, due to the faet that yi,o(i) —st xi(t), we obtam that xi(t) <st yi.n{t)- In 
the case where G{x) belongs to the class NWU, the opposite inequality holds. 

The arguments of the proof given for m = 2 remain correct for an arbitrary 
m > 2. The proof given by induction uses decomposition of the original system 
into the main system and a tagged server station as above. The further arguments 
for stochastic comparison of Markov processes are also easily extended for the case 
of an arbitrary m > 2. □ 

From the above results for the Markov processes the statement of the lemma 
follows. The stochastic inequalities between T,„ „(to — 1) and o(m — 1) follow 
by the coupling arguments. The lemma is completely proved. □ 



The results obtained in the previous section enable us to establish theorems for 
the number of losses in M/GI/m/n queueing systems during their busy periods. 

Theorem 4.1. Under the assumption A — mfj,, the expected number of losses during 
a busy period of the M/GI/m/n queueing system is the same for all n > 1. 

Proof. Consider the system M/GI/m/n under the assumption A = m/i, and sim- 
ilarly to the construction in the proof of Lemma 13.11 let us delete all the intervals 
where the number of customers in the system is less than m, and merge the corre- 
sponding ends. The process obtained is denoted ym,n(t). This is the 1-stationary 
process of orbital busy periods. 

The stationary departure process, together with the arrival 1-Poisson process of 
rate A and the number of waiting places n describes the stationary M/G/l/n queue- 
length process (with generally dependent service times). As soon as a busy period 
is finished (in our case it is an orbital busy period, see Section [2] for the definition), 
the system immediately starts a new busy period by attaching a new customer 
into the system. This unusual situation arises because of the construction of the 
process. There are no idle periods, and servers all are continuously busy. Thus, the 
busy period, which is considered here, is one of the busy periods attached one after 
another. 

Let T be a large period of time, and during that time there are K{T) busy 
periods of the M/G/l/n queueing system (which does not contain idle times as 
mentioned). Let L{T) and v{T) denote the number of lost and served customers 
during time T. We have the formula: 



the proof of which is given below. 

Relationship (|4.1|) has the following explanation. The left-hand side term EL(r)-|- 
Et^(r) is the expectation of the number of lost customers plus the expectation 
of the number of served customers during time T, and the right-hand side term 
AT -|- Eif (T) is the expectation of the number of arrivals during time T plus the 
expected number of attached customers. 

Relationship (j4.ip can be proved by renewal arguments as follows. 

There are m independent copies xj"^^(i), yi{\t), x^'"^(i) of the station- 
ary renewal process, which model the process ym^nit)- (In fact, we have m 1- 
stationary processes, which have the same distributions as m independent renewal 
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processes x'j^\t), xp'(t), . . . , x^™^(t).) Let I < i < m, and let Ci, C2,. . .CKi(T) 

be such points of busy period starts associated with the renewal process yii\t) 
(one of those m independent and identically distributed renewal processes), where 
Ki{T) denotes the total number of these regeneration point indexed by i. De- 
note also by zi, Z2, ■ • ■ , zk^^t) the corresponding lengths of busy periods, by £1, 
^2,- ■ • /Ki{T) the corresponding number of losses during these Ki busy periods, and 
by ni, 71.2,. . . ,nii.(^x) the corresponding number of served customers during these 
busy periods. Let Ti — zi + Z2 + ■ ■ ■ + ZKi{T) , let Li{T) = £1 + 12 + ■ ■ ■ + £Ki(T) and 
let yi{T) = rii + ^2 + . . . + riK^iT)- 

Since at the moments Ci, C2,. . . CKi(T) of the busy period starts the distribution 
of the above stationary Markov process of residual times is the same, then the num- 
bers of losses li, £2,- ■ ■ ,£Ki{T) and, respectively, the numbers of served customers 
ni, n2,. . . ,ri;f.(7-') during each of these busy periods have the same distributions, 
and one can apply the renewal reward theorem. 

By the renewal reward theorem we have: 

Taking into account that 

lim — — ^— ^ = m, 
lim - ^ 



and 



T^oo ELi{T) 
lim V4 1, 



Imi = 1, 



because of the correspondence between the left- and right-hand sides, from ()4.2p 
we arrive at (|4.ip . Thus, we bypass the fact that the times between departures 
are dependent, and thus (|4.ip is actually obtained by application of the renewal 
reward theorem by a usual manner, like in the case of independent times between 
departures (see e.g. Ross [28^, Karhn and Taylor [H]). 
Together with (|4.ip we have 

(4.3) lim — — — — Ez^(r) = lim — — ——miiT. 

Let us now introduce the following notation. Let C„ denote the length of an 
orbital busy period, and let L„ and correspondingly denote the number of lost 
and served customers during that orbital busy period. Using the arguments of [5], 
we prove that EL„ = 1 for all n > 1. 

Indeed, from (|4.ip and (|4.3I) we have the equations: 

(4.4) EL„ + Em„ = AEC„ + 1, 

(4.5) Em„ = m^EC„. 

The substitution A — mfj. into the system of equations (|4.4p and (|4.5p yields EL„ = 
1. 
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Hence, during an orbital busy period there is exactly one lost customer in average 
for any n > 0. To finish the proof we need in a deeper analysis. First, we should 
find the expected number of queueing periods during one orbital busy period. For 
this purpose, one can use the similar construction by deleting all the intervals 
when the number of customers in the system is not greater than m, and merge the 
corresponding ends. The obtained process is denoted ym+i,n(i), and this is one 
stationary process of queueing periods following one after another. 

The structure of the process ym+i,ri(^) is similar to that of the process Ym.n{i)- 
The process ym+i,n{t) describes the stationary M /G/l/n — 1 queueing system, the 
service times of which are generally dependent. As soon as one busy period in 
this system is finished, a new customer starting a new busy period is immediately 
attached into the system. Thus, the only difference between the processes Ym.nit) 
and ym+i.n{t) is that the numbers of waiting places differ by value of parameter 
n. Therefore, using the similar notation and arguments, one arrive at the conclu- 
sion that the expected number of losses per queueing period is equal to 1 as well. 
Therefore, in long-run period of time, the number of queueing periods and orbital 
busy periods is the same in average. So, there is exactly one queueing period per 
orbital busy period in average. 

Therefore, during a long-run period the number of events that the different 
Markov processes ym-i,n{t) change their values after deleting queueing periods 
and merging the ends (as exactly explained in the proof of Lemma [ST} is the same 
in average for all rt > 1, and the stationary characteristics of all of these Markov 
processes ym-i.n{t), given for different values n=l,2,..., are the same. Hence, 
the expectation ETm,n{'m — 1) is the same for all n—1,2,. . .as well. (Recall that 
Tm,n{'ni — 1) denote the total time during a busy period when m — 1 servers are 
occupied.) 

Hence, using Wald's identity connecting ETm^„(rn — 1) with ELm,n (the ex- 
pected number of losses during a busy period) we arrive at the desired result, since 
^Tm.ni'nT' — 1) and the expectation of the number of orbital busy periods during a 
busy period of M/GI/m/n both are the same for all n > 1. □ 

Application of Lemma 13.11 and the arguments of Theorem 14. II enables us to prove 
the following result. 

Theorem 4.2. Let A — m/i. Then, under the assumption that G{x) belongs to the 
class NBU, for the number of losses in M/GI/m/n queueing systems, n > I, we 
have 

cm^ 

(4.6) EL„,„ = -, 

m! 

where the constant c > 1 depends on m and the probability distribution G{x) but is 
independent of n. 

Under the assumption that G{x) belongs to the class NWU we have (|4.6p but 
with constant c < 1 . 

Proof. Notice first, that for the expected number of losses in M/GI /m/0 queueing 
systems we have 

_ "^^ 

i^^m.O — T 

to! 

This result follows immediately from the Erlang-Sevastyanov formula [29j . so that 
the expected number of losses during a busy period of the M/GI/m/0 queueing 
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system is the same that this of the M/M/m/Q queueing system. The expected 
number of losses during a busy period of the MjM /m/Q queueing system, ELm,o — 
is also derived in Section 
In the case where G{x) belongs to the class NBU according to Lemma 13.11 we 
have ETm^„(rn — 1) > ETm_o('7t — 1), and therefore, the expected number of orbital 
busy periods in the M/GI /m/n queueing system {n > 1) is not smaller than this 
in the M/GI /m/O queueing system. Therefore, repeating the proof of Theorem l4.1l 
leads to the inequality ELm,n > ELm,o and consequently to the desired result. If 
G{x) belongs to the class NWU, then we have the opposite inequalities, and finally 
the corresponding result stated in the formulation of the theorem. □ 

5. Batch arrivals 

The case of batch arrivals is completely analogous to the case of ordinary (non- 
batch) arrivals. In the case of a Markovian /M/m/n queueing system one can 
also apply the level-crossing method to obtain equations analogous to (|2.6p - (|2.13p . 
The same arguments as in Sections [3] and |4] in an extended form can be used for 
/GI/m/0 queueing systems. 
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